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While bare diagrammatic series are merely Taylor expansions in powers of interaction strength, 
dressed diagrammatic series, built on fully or partially dressed lines and vertices, are usually con¬ 
structed by reordering the bare diagrams, which is an a priori unjustified manipulation, and can even 
lead to convergence to an unphysical result [Kozik, Ferrero, and Georges, PRL 114, 156402 (2015)]. 
Here we show that for a broad class of partially dressed diagrammatic schemes, there exists an 
action S ^ depending analytically on an auxiliary complex parameter £, such that the Taylor ex¬ 
pansion in £ of correlation functions reproduces the original diagrammatic series. The resulting 
applicability conditions are similar to the bare case. For fully dressed skeleton diagrammatics, ana- 
lyticity of is not granted, and we formulate a sufficient condition for converging to the correct 
result. 

PACS numbers: 71.10.Fd, 02.70.Ss 


Much of theoretical physics is formulated in the lan¬ 
guage of Feynman diagrams, in various fields such as 
condensed matter, nuclear physics, and QCD. A pow¬ 
erful feature of the diagrammatic technique, used in each 
of the above fields, is the possibility to build diagrams 
on partially or fully dressed propagators or vertices, see, 
e.g., Refs. [lj-Q. In quantum many-body physics, no¬ 
table examples include dilute gases, whose description 
is radically improved if ladder diagrams are summed up 
so that the expansion is done in terms of the scattering 
amplitude instead of the bare interaction potential, and 
Coulomb interactions, which one has to screen to have a 
meaningful diagrammatic technique. 

With the development of Diagrammatic Monte Carlo, 
it becomes possible to compute Feynman diagrammatic 
expansions to high order for fermionic strongly corre¬ 
lated quantum many-body problems M- The number 
of diagrams grows factorially with the order, even for a 
fully irreducible skeleton scheme [12]. Nevertheless, for 
fermionic systems on a lattice at finite temperature, dia¬ 
grammatic series (of the form a n with a n the sum of 
all order-n diagrams) are typically convergent in a broad 
range of parameters, due to a nearly perfect cancellation 
of contributions of different sign within each order, as 
proven mathematically [131 ] and seen numerically 


ill- 


One might think that partial or full renormalization of 
diagrammatic elements (propagators, interactions, ver¬ 
tices, etc.) always leads to more compact and better be¬ 
having diagrammatic expansion. However, such a dressed 
diagrammatic series cannot be used blindly: Even when 
it converges, the result is not guaranteed to be correct, 
since it is a priori not allowed to reorder the terms of a 


series that is not absolutely convergent (the sum of the 
absolute values of individual diagrams is typically infi¬ 
nite, due to factorial scaling of the number of diagrams 
with the order). And indeed, for a skeleton series, i.e., a 
series built on the fully dressed propagator, convergence 
to a wrong result does occur in the case of the Hub¬ 
bard model in the strongly correlated regime near half 
filling 14], and preliminary results suggest that the cor¬ 
responding self-consistent skeleton scheme converges to a 
wrong result as a function of the maximal self-energy di¬ 
agram order J\T [l5j. Both of these phenomena are clearly 
seen in the exactly solvable zero space-time dimensional 

case [ 11113 - 


In this work, we establish a condition that is neces¬ 
sarily violated in the event of convergence to a wrong 
result of the self-consistent skeleton scheme. Further¬ 
more, we show that this convergence issue is absent for 
a broad class of partially dressed schemes. In particular, 
we propose a simple scheme based on the truncated skele¬ 
ton series. The underlying idea is to construct an action 
S^' 1 that depends on an auxiliary complex parameter £ 
such that the Taylor series in £ of correlation functions 
reproduces the dressed diagrammatic series built on a 
given partially or fully dressed propagator. This makes 
the dressed scheme as mathematically justified as a bare 
scheme, provided S ® is analytic with respect to £ and 
S^ =lS> coincides with the physical action; these condi¬ 
tions hold automatically in the partially dressed case, 
while in the fully dressed case they hold under a simple 
sufficient condition which we provide. Our construction 
applies to a general class of diagrammatic schemes built 
on dressed lines and vertices, including two-particle lad¬ 
ders and screened long-ranged potentials. 
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Partially dressed single-particle propagator. We con¬ 
sider a generic fermionic many-body problem described 
by an action 

S[tp,ljj} = (V’lG'oV) -t-S'int bP,$\ (1) 

where if .!>, if) are Grassmann fields 0 , and we use bra- 
ket notations to suppress space, imaginary time, possi¬ 
ble internal quantum numbers, and integrals/sums over 
them, i.e., (V’|Gq 1 |V’) denotes the integral/sum over r, 
r and o of ^(r, t) (Gg \ ^ <T )(r, r). Gg 1 stands for the 
inverse, in the sense of operators, of the free propaga¬ 
tor. The full propagator G and the self-energy E are 
related through the Dyson equation G -1 = Gg 1 — E. 
The bare Feynman diagrammatic expansion corresponds 
to perturbation theory in Si nt . In order to generate a di¬ 
agrammatic expansion built on a partially dressed single¬ 
particle propagator Gw, we introduce an auxiliary action 
of the form 

SjpbP,4>\ = (^|Ggjv(£)| V>) +^S int [il),if)], (2) 

where 

^o.aKO = + £Ai + • • • + C^Aw, (3) 

£ is an auxiliary complex parameter, and Ai,..., Aw are 
appropriate operators. Gw is the single particle prop¬ 
agator for ^- 0) . At £ ^ 0, one can still view Gy as 
the free propagator, provided one includes in the interac¬ 
tion terms not only fS, n t, but also the quadratic terms 
(ip\t; n A n \ip). Accordingly, £ is interpreted as a coupling 
constant, and the £"A ra acquire the meaning of counter¬ 
terms. These counter-terms can be tuned to cancel out 
reducible diagrams, thereby enforcing the dressed char¬ 
acter of the diagrammatic expansion. A natural require¬ 
ment is that S%- 1] coincides with the physical action S , 
i.e., that 

AT 

Gjy 1 + X] A" = G 0 1 . (4) 

n =1 

For given Go, this should be viewed as an equation to 
be solved for Gy (it is non-linear if the A„’s depend on 
Gw). The unperturbed action for the dressed expansion, 
(^IGj^ | if), is shifted by the A„’s with respect to the 
unperturbed action for the bare expansion, (iflG^lif). 

We can then use any action of the generic class (0 for 
producing physical answers in the form of Taylor expan¬ 
sion in powers of £, provided the propagator Gw and the 
shifts A„ satisfy Eq. 0 . More precisely, consider the full 
single-particle propagator Gw(£) °f the action Sffi, and 
the corresponding self-energy 

^(e):=G 0 -^(0-G^(0. (5) 

Note that since £w _1) = we have Gw(£=l) = G and 
hence also Ew(£=l) = 5k We assume for simplicity that 


Ew(£) is analytic at £ = 0, and that its Taylor series 
^w^ [Gw] £ n , converges at £ = 1. We expect these 
assumptions to hold for fermionic lattice models at finite 
temperature in a broad parameter regime, given that the 
action sffi is analytic in £ Then, 

since S$ is an entire function of £, we can conclude that 

OO 

E = £e$[Gw], (6) 

n= 1 

i.e., the physical self-energy is equal to the dressed dia¬ 
grammatic series. 

This last step of the reasoning can be justified using 
the following presumption: Let P be a connected open 
region of the complex plane containing 0. Assume that 
is analytic in V, that the corresponding self-energy 
E(£) is analytic at £ = 0, and that E(£) admits an an¬ 
alytic continuation E(£) in V. Then, E and E coincide 
on V. This presumption is based on the following ar¬ 
gument: Since S^ is analytical, if no phase transition 
occurs when varying £ in T>, then E(£) is analytical on T>, 
and by the identity theorem for analytic functions, E and 
E coincide on V. If a phase transition would be crossed 
as a function of £ in T>, analytic continuation through the 
phase transition would not be possible |20], contradicting 
the above assumption on the existence of E. Applying 
this presumption to E(£) := 5qy^ [Gw] £ n , which 

has a radius of congergence R > 1 (from the Cauchy- 
Hadamard theorem), and taking for V the open disc of 
radius R, we directly obtain Eq. 0 provided R > 1. If 
R= 1, we can still derive Eq. 0, using Abel’s theorem 
and assuming that Ew(£) is continuous at £=1, which, 
given that the action is entire in £, is generically expected 
(except for physical parameters fined-tuned precisely to 
a first-order phase transition, where E is not uniquely 
defined). 

Semi-bold scheme. We first focus on the choice 

An = S^idlGw] (l<n<W), (7) 

where E^"^ [Q] is the sum of all skeleton diagrams of or¬ 
der n, built with the propagator Q and the bare in¬ 
teraction vertex corresponding to S\ nt , that remain con¬ 
nected when cutting two Q lines. This means that Gw 
is the solution of the bold scheme for maximal order Af, 
cf. Eq. 0. For a given Af, higher-order dressed graphs 
can then be built on Gw- The numerical protocol corre¬ 
sponding to this ‘semi-bold’ scheme consists of two inde¬ 
pendent parts: Part I is the Bold Diagrammatic Monte 
Carlo simulation of the truncated order- Af skeleton sum 
employed to solve iteratively for Gw satisfying Eqs. 00; 
Part II is the diagrammatic Monte Carlo simulation of 
higher-order terms, T,ff[Gjy],n > Af, that uses Gw as 
the bare propagator. Note that here Af is fixed (con- 
trarily to the conventional skeleton scheme discussed be¬ 
low) , and the infinite-order extrapolation is done only in 
Part II. 
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The Feynman rules for this scheme are as follows: 

[G^] = EWj [G M ] for n < Af-, (8) 

while for n > J\T + 1, eW[(5/v] is the sum of all bare 
diagrams, built with Gjy as free propagator and the bare 
interaction vertex corresponding to S'int, which do not 
contain any insertion of a subdiagram contributing to 
E^JG*] with n < AT. Indeed, each such insertion is 
exactly compensated by the corresponding counterterm. 
To derive Eq. (j8]), we will use the relation 

OO 

s*-(o = £ E (9) 

n— 1 


lEP(P> r )l < C 1 (p,r)- and 

(ii) Ga/(p,t) is uniformly bounded, i.e., there exists a 


constant C 2 such that for all (Af, p, r), |Gw(p, 
Our derivation is based on the action 

r)| < G 2 . 


:= EEo^ } - 

(11) 

Clearly, 

S^ 1 = (V’l G _1 + \if) + £S' int 

(12) 

with 

L (?) (p,r) := lim L$( p,r). 

AI —>00 

(13) 


where = stands for equality in the sense of formal power 
series in £, and we will show the proposition 

k 

= E E boi d [^] r+o(e fe+1 ) iv k ) 


for any k £ {0,. .., AT + 1}, by recursion over k. (Pk= o) 
clearly holds. If (Tk) holds for some k < Af , then we 
have G^(£) = G+ 0(£ fc+1 ), as follows from Eqs. (0, 
0 and (0. Substitution into Eq. 0 then yields (Vk+ 1 ). 

Alternatively to the semi-bold scheme Eq. 0, other 
choices are possible for the shifts Ai,..., A jy and the 
dressed propagator Gjy. For example, the shifts can be 
based on diagrams containing the original bare propaga¬ 
tor Go instead of Gjy. In the absence of exact cancella¬ 
tion, all diagrams should be simulated in Part II of the 
numerical protocol, and A n will enter the theory explic¬ 
itly. This flexibility of choosing the form of A„’s, along 
with the obvious option of exploring different Af’s, pro¬ 
vides a tool for controlling systematic errors coming from 
truncation of the ^-series [23[. 

Skeleton scheme. We turn to the conventional scheme 
in which diagrams are built on the fully dressed single¬ 
particle propagator. The corresponding numerical pro¬ 
tocol is identical to Part I of the above one, with the 
additional step of extrapolating Af to infinity, as done 
in [1- 11, 211. Accordingly, we assume that the ‘skele¬ 
ton sequence’ G^f converges to a limit G when Af —> oo. 
The crucial question is under what conditions one can 
be confident that G is the genuine propagator G of the 
original model. The answer comes from the properties of 
the sequence of functions 


4 ): =E41[w 


( 10 ) 


Let us show that G = G holds under the following suffi¬ 
cient condition: 

(i) for any £ in a disc V = {|£| < R} of radius R > 1, 
and for all (p, r), L^(p,r) converges for Af^-oo; more¬ 
over this sequence is uniformly bounded, i.e., there ex¬ 
ists a function Gi(p,r) such that V£ € V, V(A/", p, r), 


Since Sjy- ^ = S, we have Sx ^ = S, and thus 
Goo(£=l) = G where G«,(£) is the full propagator of 
the action S& . 

We first observe that (p, r) is an analytic function 
of £ £ V for all (p, r), and that 


1 d n 
n\ <9£ n 


L^( p, 


T ) 


£=o 


y(«) 

^bold 


[G](p,r). 


(14) 


This follows from conditions (i.ii), given that momenta 
are bounded for lattice models. Indeed, for any trian¬ 
gle T included in V , § r d£ L$ (p, r) = 0. Thanks to 
condition (i), the dominated convergence theorem is ap¬ 
plicable, yielding § r d£ L^> (p, r) = 0. The analyticity of 
£ i—^ follows by Morera’s theorem. To derive 

Eq. (fTdl) we start from 


I EL rl 

n\ « 9 £» L Af(P> ) 




(15) 


S=o 


By Cauchy’s integral formula, the l.h.s. of Eq. m 
equals 1/(2*7 t) <f c dfL$(p, r)/£ n+1 where C is the unit 
circle. Using again condition (i) and the domi¬ 
nated convergence theorem, when Af—too, this tends 
to 1/(2*7 t) j c dfL^(p, r)/£" +1 , which equals the l.h.s. 

of Eq. (HU). To show that Eoidl^A/’KP’ r ) tends to 
E boid[^](P, r )> we consider each Feynman diagram sepa¬ 
rately; the dominated convergence theorem is applicable 
thanks to condition (ii), the boundedness of the integra¬ 
tion domain for internal momenta and imaginary times, 
and assuming that interactions decay sufficiently quickly 
at large distances for the bare interaction vertex to be 
bounded in momentum representation. 

Hence 


L (0 = £ £("J d [G] £". (16) 

n= 1 

As a consequence, the action S E generates the fully 
dressed skeleton series built on G, i.e., its self-energy 
£oo(£) has the Taylor expansion X))LLi E b"id [^] i an d 
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the Taylor series of Goo(£) reduces to the ^-independent 
term G. This can be derived in the same way as Eq. ©> 
by showing by recursion over k that for any k > 0, 
Eoo(0 = £*=i£&] d [G]e + 0(£ fe+1 ). Furthermore, 
having shown above the analiticity of L^\ i.e., of , 
we again expect that Goo(£) is analytic at £ = 0 (for 
fermions on a lattice at finite temperature), and we 
can use again the above presumption to conclude that 
Goo(£ = 1) = G = G. 

Dressed pair propagator. So far we have discussed 
dressing of the single-particle propagator while keeping 
the bare interaction vertices. We turn to diagrammatic 
schemes built on dressed pair propagators. We restrict 
to spin- 1/2 fermions with on-site interaction: 

^ w3 

= U ^2 dr (V'tV’lV’4-V’t)( r -' r )> ( 17 ) 

r 

where U is the bare interaction strength. For simplicity 
we discuss dressing of the pair propagator while keep¬ 
ing the bare Go- It is necessary to perform a Hubbard- 
Stratonovich transformation in order to construct the ap¬ 
propriate auxiliary action. Introducing a complex scalar 
Hubbard-Stratonovich field rj leads to the action 

S[ip, tp, V, V\ = WGfW - - (r?|n 0 |r?) 

+ (! 8 ) 

where n 0 (r,r) = — (Go,fGo,4.)(r, r) and r 0 is the sum of 
the ladder diagrams, Tg 1 (p, 0„) = f7 _1 — n 0 (p, fl n ) with 
il n the bosonic Matsubara frequencies. 

We first consider the diagrammatic scheme built on Go 
and r 0 . We denote by £|”j[Go,r 0 ] the sum of all self¬ 
energy diagrams of order n, i.e. containing n To-lines. 
This diagrammatic series is generated by the shifted ac¬ 
tion 

5 iad l^,4’, r hV\ = (4 ’\Gq 1 \4>) - ~^ 2 (ri\Il 0 \r]) 

+ £ (MVm.V’t) + (V , 4-V’tl ? ?)) - ( 19 ) 

in the sense that self-energy Ei a d(£) corresponding to this 
action has the Taylor series Fijad Kug Fo] £ 2 ”. In¬ 

deed, the counter-term £ 2 n 0 cancels out the reducible 
diagrams containing GqGq bubbles. Therefore, if this 
diagrammatic series converges, then it yields the physi¬ 
cal self-energy. This follows from the same reasoning as 
below Eq. ©• The same applies to the series for the pair 
self-energy II in terms of [Go,To]. Here n is defined by 
T -1 = Tf / 1 — n, where T denotes the fully dressed pair 
propagator, used in UM- 

More complex schemes, built on other dressed pair 
propagators than To, can be generated by the shifted 
action 

+ £ (MVm.V’t) + , ( 20 ) 


where 

r o>(0 = f ^• 1 +£ 2 n 1 + ... + £ 2A ^ A r (21) 

and one imposes ro,Af(£ = 1) = To- In particular, the 
semi-bold scheme is defined by 

( 22 ) 

where ^"^[ 7 ] is the sum of all skeleton diagrams of 
order n built with the pair-propagator 7 that remain 
connected when cutting two 7 -lines. As usual, H^^ = 
—GG + Go Go. This scheme was introduced previously 
for A/"=l @. 

Finally, we consider the skeleton scheme built on Go 
and r. Assuming that the skeleton sequence f/y con¬ 
verges to some r, one can show analogously to the above 
reasoning that T is equal to the exact T under the fol¬ 
lowing sufficient condition: 

(i) for any £ in a disc V = (|£| < A} of radius R > 1, and 

for all (p, O 0 , Mjp (p, n„) := Eli I^Mp, O.) C 
converges for TV—>-oo; moreover this sequence is uni¬ 
formly bounded, i.e., there exists G(p, fl n ) such that 
€ £>, V(A/’,p,fi n ), \Mffi(p,n n )\ < C(p,fl„); and 

(ii) Tjg-(p,rt n ) is uniformly bounded. 

Screened interaction potential. Finally, we briefly ad¬ 
dress the procedure of dressing the interaction line, which 
is particularly important for long-range interaction po¬ 
tentials. Restricting for simplicity to a spin-independent 
interaction potential V(r), the interaction part of the ac¬ 
tion writes 

7^51 J dT - r’) ,t). 

r.r' ^ 

(23) 

We again keep the bare Go for simplicity and consider 
dressing of V only. Introducing a real scalar Hubbard- 
Stratonovich field x leads to the action 

s [= (V , |GoV) + ^(x|i /_1 lx) + * J2{x\^a)- 

(24) 

Here we assume that the Fourier transform V (q) of the 
interaction potential is positive, so that the quadratic 
form (x| V^ 1 |x) = (27r)- d Jjf dr J d d q |%(q, r)| 2 /H(q) is 
positive definite. The auxiliary action takes the form 

Sy x] = (i/IGq 1 |V ; ) 

+ 2 (x|Fa/ 1 +C 2 H i + ■ • .+£ 2 A ^Hy|x) + *£ 

(25) 

The semi-bold scheme corresponds to = ll[”] d [ Vy ] 
where n now stands for the polarization. In particular, 
V\ is the RPA screened interaction. 
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Summarizing, we have revealed an analytic structure 
behind dressed-line diagrammatics. More precisely, we 
have exhibited the function which analytically continues 
a dressed diagrammatic series. This function originates 
from an action that depends on an auxiliary parame¬ 
ter When the action is a polynomial in £, the situ¬ 
ation reduces to the one of a bare expansion. Within 
this category, a particular case well suited for numerical 
implementation is the senri-bold scheme for which the 
bare propagator is taken from the truncated bold self- 
consistent equation. For the fully bold scheme, we con¬ 
struct an appropriate auxiliary action, but only under 
a certain condition. If this condition is verified numeri¬ 
cally, it is safe to use the fully bold scheme. If not, the 
semi-bold scheme remains applicable. 

Furthermore we have demonstrated the generality of 
the shifted-action construction by treating the case of 
a dressed pair propagator and of a screened long-range 
interaction. Further extensions left for future work are 
dressing of three-point vertices, as well as justifying re¬ 
summation of divergent diagrammatic series by consid¬ 
ering non disc-shaped analyticity domains D. 
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